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Abstract 

The Chaplygin gas model, characterized by an equation of state of the 
type p = emerges naturally from the Nambu-Goto action of string 
theory. This fluid representation can be recast under the form of a tachy- 
onic field given by a Born-Infeld type Lagrangian. At the same time, the 
Chaplygin gas equation of state can be obtained from a self-interacting 
scalar field. We show that, from the point of view of the supernova type la 
data, the three representations (fluid, tachyonic, scalar field) lead to the 
same results. However, concerning the matter power spectra, while the 
fluid and tachyonic descriptions lead to exactly the same results, the self- 
interacting scalar fleld representation implies different statistical estima- 
tions for the parameters. In particular, the estimation for the dark matter 
density parameter in the fluid representation favors a universe dominated 
almost completely by dark matter, while in the self-interacting scalar field 
representation the prediction is very closed to that obtained in the ACDM 
model. 
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1 Introduction 



The Claplygin gas equation of state, given by 



was first presented in the reference [T] in the context of the study of incompress- 
ible fluids. Later, it has revealed important in the study of the aerodynamics 
problem [2!|. More recently, it has re-appeared as an example of a irrotational 
fluid that can be supersymmetrized [3]. Moreover, it has been discovered an 
interesting connection between this equation of state and the Nambu-Goto ac- 
tion of string theory written in the light-cone variables [4]. Its application to 
cosmology was first pointed out in [5J. A phenomenological generalization has 
been proposed [6] , such that the equation of state takes the form 



where a is a free parameter, generally a positive number. 

The interest for the Chaplygin gas model, and its phenomenological gener- 
alization, in cosmology emerged from the evidences in favor of an accelerating 
expansion of the universe at present time [3 [8] . These evidences ask for a fluid 
with negative pressure that begins to dominate the matter content of the uni- 
verse at about z ^ 1 {z being the cosmological redshift), whithout affecting 
the previous history of the universe. In particular the primordial nucleosynthe- 
sis, the transition from a radiation dominated universe to a matter dominated 
universe and, finally, the process of structure formation driven by cold dark 
matter must remain unaltered. In spite of being a natural candidate for the 
fluid responsible for the cosmic acceleration, the cosmological constant faces 
many theoretical and observational problems, mainly represented by the huge 
discrepancy of the observed and the predicted values, as well as the necessary 
flne tunning of its value to achieve an acceleration just as measured at present 
time. 

The Chaplygin gas has revealed an interesting alternative to the so-called 
ACDM model. In fact, considering a Friedmann-Lemaitre-Robertson- Walker 
universe (FLRW) , a fluid with energy density p and pressure p must satisfy the 
conservation law 



where a is the scale of factor of the spatial section. For the equation of state 
([1]), this equation can be easily integrated leading to 



where B is an integration constant. The important point concerning gl) is the 
fact that it contains relevant asymptotic regimes: as a — > 0, i.e. as we go back to 



p + 3-{p + p)=0, 



(3) 




(4) 
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the past, p ^ a~^, reproducing a pressureless fluid (cold dark matter); moreover, 
as a — > cxD, p ^ constant, leading to a behaviour typical of a cosmological con- 
stant. Hence, the transition from a non-accelerating to an accelerating universe, 
can be achieved dynamically. At same time, the process of structure formation 
is guaranteeded by the fact that the Chaplygin gas mimics cold dark matter at 
previous time. 

The confrontation of the Chaplygin gas model with observations have shown 
that it is quite competitive with other dark matter/dark energy models, like the 
ACDM and quintessence models. However, there are some interesting aspects 
that point out to the necessity of more deep investigations. For example, while 
the confrontation with the supernova type la (SN la) data shows that the unifi- 
cation model is preferred [9] , that is essentially no exotic new component besides 
the Chaplygin gas, the structure formation analysis reveal that a large amount 
of extra dark matter is still needed [101 HI]- Moreover, supernova data seems 
to favor, in the case of the generalized Chaplygin gas model, negatives values 
of the parameter a. If we turn now structure formation, the study of the cases 
with negative values of a becomes more sensible since this implies imaginary 
sound velociy, hence plagued with instabilities [12]. Finally, mainly in the case 
of structure formation, the results may be sensitive on how the Chaplygin gas is 
represented, that is, if it is viewed as a tachyonic field, a fluid or a self-interacting 
scalar field. 

Our goal here is to consider this last problem: how the use of different repre- 
sentations of the Chaplygin gas affects the final estimation of the cosmological 
parameters. We use the SN la data and the 2dFGRS data for the matter power 
spectrum. In special, we will show that the tachyonic and fluid representations 
give the same results even for the matter power spectrum. However, the self- 
interacting scalar field, even if it admits in principle an extension for negative 
values for the parameter a when the generalized Chaplygin gas is considered, 
leads to very different parameter estimations. 

From now on, for the sake of simplicity, we will work only with the tradi- 
tional Chaplygin gas model defined by the equation of state H]). The tachyonic 
formulation can be obtained from the Lagrangian 



where T is the tachyonic field. The corresponding energy-momentum tensor is 
then given by 



For the case of the FLRW fllat universe, this energy-moment tensor leads to the 
following expressions for the energy density and for the pressure: 




(5) 



V[T) 



X^T, + V{T)^l-T,pT'Pg^; 



(6) 



^1 - T-pT-P 



V{T) 



(7) 



PT = 



PT 




(8) 
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Thus, equation ^ with A = constant imphes that V{T) must be a constant, 
that is, independent of T. For some interesting generahzations with a non- 
constant V{T) see, for example, [IB] . 

Let us consider now the energy-momentum tensor of a perfect fluid, 

Tfj^i^ = {p + p)ui_,Uu - pg,iu, (9) 

where is the four-velocity of the fluid. Introducing small perturbations 
around a background value (dp, Sp, 5u^ and Sg^^i, — using the synchronous 
coordinate condition h^Q — 0, we obtain at flrst order the following expressions: 

6T(3Q = 5p, STi-j ^ -phij - Spg^j . (10) 

Using now the expression ([5]), we obtain 

VtSt 

oTm = =^ opT, (11) 



(1 - r2)3/: 

vfdf 



STij = , . ffij + vV 1 - t'^hij = -prhij - Sprgij- (12) 

VI -T2 

Hence, the fluid approach represented by ^ reproduces the features of the 
tachyonic model described by (O at the background and at the perturbative 
levels. 

Since the tachyonic fluid is entirely acquainted by the fluid represntation 
(at classical level, of course), let us consider the Einstein's equation coupled to 
baryonic matter and the Chaplygin gas. The fleld equations read, 

R^,=87rG(^T'^,-^g^,T''^ + SvrG (^T^, - ig,,,T^^ , (13) 

rr,, = 0, rr;^ = 0. (14) 

The superscript (subscript) b and c stands for baryons and Chaplygin. 

We perturb now these fleld equations, using again the synchronous coordi- 
nate conditon. We divide all equations by the Hubble parameter today Hq in 
order to make the equations dimensionless. We end up with the following set of 
coupled equations for the perturbed quantites: 

■i o^r 3^ ^ 3^ 4ia6 + l~A^ , . 

djn + - -i'lmOm = 7:"c T r , tTC'c, (15) 

a 2 2 Aa° + [1 — A) 

f 1-^ Aa^ , 1-^ . 

Aa^ + l-A '^-^'a Aa^ + l-A ^^ ^ Aa^ + 1 - A ^^ (^'^ 

0.+ ^- ^i:'^'-f 0. ^ ^kH.. (17) 
a Aa^ + 1- A I- A ^ ' 

In these expressions 5ra and 5c are the density contrast for baryons and for 
the Chaplygin gas respectively, while 9c is the perturbation in the four velocity 
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Fi glire 1 '. Two-dimensional probability distribution for and A using only the 2dFGRS 
data (left), only the SN la data (center) and the joint probability from both sets of data. The 
lighter is the region, the higher is the associated probability density. 



connected with the Chaplygm gas. We have noted A — -4-, which is the sound 

PaO 

velocity of the Chaplygin gas fluid, with pco = ^A + B being the corresponding 
energy density today. 

Passing now from the (dimensionless) time variable t to the scale factor 
variable a, the above equations take the following form: 

,, l-Ae^ Aa^ 1 - A ^ ^ 

Aa^ + l-Aj + '^Aa^ + l-A^^ = Aa^ + l^A^-^ ^^^^ 
^, , Aa^ + 2{l-A) 9, Aa^ kHl 

Aa^ + l-A 7 ^ ^^^^ 

The primes mean derivative with respect to a and Iq is the Hubble radius today. 
Moreover, we have the following definitions: 

/(a) = + n,{a)a, (21) 

a„(a) n,{a)l-A-2Aa'' 
q(a) — — a 1^^) (22) 



= n,o]lA + ^—^. (24) 

Note that the sound speed of the fluid plays an important role, in particular 
multiplying the term originated from the gradient of pressure in equation (j20p . 
It is this fact that leads to instabilities problems when the generalized Chaplygin 
gas model, with a negative, is considered. 
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The analysis for the SN la data is made through the computation of the 
luminosity distance, given by 



c dz' 
Dl^-^{1 + z) ^ (25) 

^^o' 'Jo V^WITTW) 

where z is the redshift given by z = — 1 + ^ . 

We note that, concerning the matter power spectrum we have two free pa- 
rameters: flmo (or equivalently QcOi since flmo + ^co = 1) and A. For the SN 
la analysis we have one more free parameter: the Hubble reduced parameter 
h defined by Hq = 100 h km/ Mpc.s. We use the statistical function which 
gives the quality of the fitting of the observational data by the theoretical model 
(see [9]). Marginalizing the probability distribution function for the SN la data 
in the variable h (it means, integrating in this variable), we end up with two 
sets of observational constraints, for the 2dFGRS and SN la data, depending 
on the parameters A and r^mo- 

We compute now the probability distribution function using the 2dFGRS 
power spectrum data. The initial conditions are fixed using the prescription 
described in reference [14j . The luminosity distance is also computed, for the 
gold sample, by using the method sketched in 9 . The results are shown in 
figure 1, where the two-dimensional PDF using only power spectrum data, only 
the SN la data and the joint probability are displayed. In figure 2, the corre- 
sponding estimations for the one-dimensional PDF for the matter density today 
are shown, while in the figure 3 the one-dimensional PDF for the sound velocity 
A for the three sets of data is displayed. The results for the SN la data agree 
with those described in reference [9]: the unified scenario, with essentially no 
dark matter, is preferred; at the same time, very large values for A, that is, 
the " quasi" -ACDM models, are also preferred. Using the 2dFGRS data, the 
conclusions of references [TOl [11] are re-obtained: the anti-unified model, with 
essentially all matter in the form of dark matter, is preferred; for the parameter 
A, again large values are preferred, but there is also a high probability distribu- 
tion near zero. The combination of both data privilegizc the scenarios indicated 
by the 2dFGRS even if some new effects appear like the maximum value for 
the PDF of no for the joint SN la and 2dFGRS data in figure 2. But there 
is some subtle points concerning this: if the joint higher dimensional PDF in 
the parameter space is first constructed, as we have done in the present work, 
than the results arc those displyaed in the figure 3; however, if the correspond- 
ing one-dimensional PDF's are combined after the marginalization on the other 
variables, the conclusions are very different, in the sense that the SN la data 
are favored. But, we think the first procedure is more consistent. 



Now, we look for a self-interacting scalar field that mimics the Chaplygin 
gas model in presence of a baryonic component. The equations of motion are 
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Figure 3: One-dimensional probability distribution for A using only the 2dFGRS data (left), 
only the SN la data (center) and the joint probability from both sets of data. 



now the following: 



„a fa 
2-+ - 

a \a 



a 



(26) 

(27) 
(28) 



This system of equations is equivalent to the system (|13I14|) if the derivative of 
the (j) and the potential V are given in terms of the scale factor as 



A 



l-A 



1/4 ' 



V 



A- 



1~A 



(29) 



(30) 



These relations allow to define implicitly V as function of cj). In fact, using the 
Friedmann equation, it is possible to integrate (|29p. However, apparently there 
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is no way to have a closed expression for V as function of 0, except in the case 
the baryonic component is absent [5lll5|. But, for our purposes, all we need are 
the relations (|29l30p . 

We perturb now the Einstein-scalar field system using again the synchronous 
gauge condition. Following the same steps as before, and re-expressing the 
final equations in terms of the variable a (the scale factor) we end up with the 
following set of coupled equations: 

In these expressions, $7™, / and g are defined as before, l{a) ~ cj), with (/) given 
by and A = Scj). Moreover, we have the following definitions: 

V4, = ^K, (33) 

f2 f ( f I'f] 

V^4' = 7J^- + f|^--^jK. (34) 

The subscript a indicates derivatives with respect to the variable a. 

This system has the advantage of admitting a straighforward extension to 
the generalized Chaplygin gas model without present instability problems when 
a is negative. In fact, as in the case of any usual scalar theory, the Laplacian 
term that generates the dependence on the wavenumber k does not contain the 
sound velocity, as it happens in the fluid reprsentation. However, at the same 
time, this scalar version of the Chaplygin gas model leads to complete different 
predictions concerning the parameter estimations. The figures 3, 4 and 5 show 
the two-dimensional and onc-dimcnsional PDF for the parameters ilmo and A, 
using only the 2dFGRS or the SN la data, and when the joint probability are 
considered. The SN la results are the same as before, since it is sensitive to the 
background only. But, the results for the matter parameter density using the 
2dFGRS are completely different: they predict a small value for flmo, as it is 
obtained for the ACDM or quintessence models, together with a smaller pic of 
probability near flmo = 1- The 2dFGRS data favors a small, but non null, value 
for A while the SN la favors a value of A near 1. The joint probability gives a 
maximum value for the probability for values of A near 0.8. 



Hence, the scalar version of the Chaplygin gas model keeps the property of 
the original formulation, in terms of the tachyonic field, only for the background. 
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Figure 4: Two-dimensional probability distribution for flo and A using only the 2dFGRS 
data (left), only the SN la data (center) and the joint probability from both sets of data. The 
lighter is the region, the higher is the associated probability density. 




Figure 5: One-dimensional probability distribution for f!o using only the 2dFGRS data 
(loft), only the SN la data (center) and the joint probability from both sets of data. 
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Fi gure 6; One-dimensional probability distribution for A using only the 2dFGRS data (left), 
only the SN la data (center) and the joint probability from both sets of data. 

As a consequence it seems that in the analysis of the structure formation, and 
considering the generahzed Chaplygin gas model, we must consider a space 
parameter where only the positive values of a are allowed. This may force 
also to restrict the SN la analysis to this region of the space parameter. Of 
course, the final estimations for the cosmological parameter sare modified by 
this restriction (see discussion in [3). We hope to present a full analysis of this 
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issue in a future paper. 
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